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O : ABSTRACT 

o 

o I We report the discovery of the two-image gravitationally lensed quasar SDSS 

J133222.62+034739.9 (SDSS J1332+0347) with an image separation of A9 = 
1'!14. This system consists of a source quasar at z s = 1.445 and a lens galaxy at 



00 



C\l I Zi — 0.191. The agreement of the luminosity, ellipticity and position angle of the 

lens galaxy with those expected from lens model confirms the lensing hypothesis. 
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1. Introduction 

Strong gravitational lensing is a powerful probe of the large scale properties of the 
universe (Kochanek 2006). The well-known physics makes it straightforward to use lensing 
in cosmological estimates of the dark energy (Turner 1990; Fukugita et al. 1990) and the 
Hubble constant (Refsdal 1964). The gravitational lensing effect depends solely on the mass 
density of the field, it is also an ideal tool for studies of dark matter in galaxies and clusters 
of galaxies (e.g., Oguri et al. 2004a). The Cosmic Lens All Sky Survey (CLASS; Myers et 
al. 2003; Browne et al. 2003), which was conducted in the radio band, is the largest current 
statistical sample of lensed compact sources. However, the CLASS has only 13 lensed radio 
sources in its well-defined statistical sample. Larger lens samples are critical to using lensed 
quasars for cosmological and astrophysical tests. 

Optical surveys for lensed quasars are both complementary to radio surveys like the 
CLASS and likely to produce larger lens samples since only one-tenth of quasars are radio- 
loud (Ivezic et al. 2002). Our Sloan Digital Sky Survey Quasar Lens Search (SQLS; Oguri et 
al. 2006) is based on the Sloan Digital Sky Survey (SDSS; York et al. 2000), and its catalog of 
roughly 100,000 spectroscopically-confirmed quasars (e.g., Schneider et al. 2005). So far, the 
SQLS has succeeded in discovering 14 new lensed quasars (Inada et al. 2003a,b,c, 2005, 2006; 
Johnston et al. 2003; Morgan et al. 2003; Pindor et al. 2004, 2006; Oguri et al. 2004b, 2005; 
Buries et al. 2006), as well as recovering many previously known lensed quasars (Walsh et al. 
1979; Weymann et al. 1980; Surdej et al. 1987; Bade et al. 1997; Oscoz et al. 1997; Schechter 
et al. 1998; Morgan et al. 2001). Indeed, the number of SDSS-discovered lensed quasars is 
already a significant fraction of all known quasar lenses (~ 80 2 ). The large number of lenses 
contained in the SQLS and its well-defined selection function means that the lens catalog 
provided by the SQLS will be very useful for statistical studies of strong lensing (Oguri et 
al. 2006). 

In this paper, we report the discovery of a gravitationally lensed quasar, SDSS J133222. 62+034739. 9 
(hereafter SDSS J1332+0347). The quasar was confirmed to be doubly imaged by a low- 
redshift early-type galaxy, from observations using the Subaru 8.2-meter telescope and the 
University of Hawaii 2.2-meter (UH88) telescope. This paper is organized as follows. The 
SDSS data for SDSS J1332+0347, including its selection, is described in §2. The follow-up 
images and spectra are shown in §3 and §4, respectively. In §5 we model the lens. Finally 
we summarize our results in §6. Throughout the paper we assume a cosmological model 
with the matter density Qm = 0.27, cosmological constant Q\ = 0.73, and Hubble constant 
h = ifo/100 kmsec _1 Mpc _1 = 0.7. Magnitudes are measured in AB system. 



2 see CASTLES website, http://cfa-www.harvard.edu/castles/ 
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2. SDSS Data 

The SDSS is a survey to conduct both optical imaging and spectroscopy. It covers about 
10,000 square degrees of the sky approximately centered on the North Galactic Cap with the 
dedicated wide-field 2.5-m telescope (Gunn et al. 2006) at the Apache Point Observatory 
in New Mexico, USA. The five broad-band photometric data are taken in imaging surveys 
(Fukugita et al. 1996; Gunn et al. 1998; Lupton et al. 1999; Tucker et al. 2006) and are reduced 
automatically by the photometric pipeline (Lupton et al. 2001). From the magnitudes, colors, 
and morphologies of each object, quasar and galaxy candidates are selected as spectroscopic 
targets (Eisenstein et al. 2001; Richards et al. 2002; Strauss et al. 2002), and then assigned for 
spectroscopic observations by the tiling algorithm (Blanton et al. 2003). The spectra cover 
the wavelength range from 3800 A to 9200 A with a resolution of R ~ 1800. The astrometric 
accuracy is better than ~ O'.'l rms (Pier et al. 2003) and the photometric zeropoint errors are 
less than ~ 0.03 magnitude over the entire survey area (Hogg et al. 2001; Smith et al. 2002; 
Ivezic et al. 2004). Most of the data have already become publicly available (Stoughton et 
al. 2002; Abazajian et al. 2003, 2004, 2005; Adelman-McCarthy et al. 2006). 

SDSS J1332+0347 was selected as a lensed quasar candidate from the SDSS quasar 
sample, using the same selection algorithm as that described in Inada et al. (2003a). In 
this algorithm, the morphological parameters of each quasar are used to select lensed quasar 
candidates. SDSS J1332+0347 was also selected by an improved algorithm (Oguri et al. 
2006) whose selection function has been extensively tested, and it is in the well-defined 
source quasar sample of the SQLS for lensed quasar statistics. 

Figure 1 shows the SDSS i-band image of the field under the seeing of 1'.'6. The SDSS 
Point Spread Function (PSF) magnitudes and the quasar redshift are summarized in Table 
1. The SDSS images indicate that the object is more extended than a single star, making it 
a good lens candidate and a target for higher angular resolution observations. 

3. Imaging Observations 

We obtained a 90 sec i-band image of SDSS J1332+0347 with the Suprime-Cam (Miyazaki 
et al. 2002) on the Subaru 8.2-meter telescope at Mauna Kea in seeing 0'.'5 on 2003 May 28. 
The i-band image was reduced in a standard way for the Suprime-Cam data using the NEKO 
software (Yagi et al. 2002) and the SDFRED package (Ouchi et al. 2004). The absolute flux 
was calibrated using nearby stars in the SDSS catalogs. The small differences between the 
SDSS and Suprime-Cam i-band response functions are not important here. We also acquired 
a 720 sec if-band image of this object with the QUick near-InfraRed Camera (QUIRC) on 
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the UH88 telescope in Cf/8 seeing on 2005 February 20. The if -band image was reduced in a 
standard way using IRAF 3 , and calibrated using the standard star FS 23 (Hawarden et al. 
2001). 

The i and H images are shown in Figure 2. Taking a closer look at the images, especially 
the i-band Suprime-Cam image, we see two compact sources bracketing an extended elliptical 
source. Indeed, models of the images using a public software, GALFIT (Peng et al. 2002), 
consisting of two point sources and an extended object with a de Vaucouleurs profile between 
them fit the data well. We used stars in the images as templates for the PSF in this procedure. 
We also show the observed galaxy after subtracting models for the two point sources, the 
observed stellar components after subtracting models for the galaxy, and the residuals after 
subtracting models for the point sources and the galaxy, in the second, third, and bottom 
panels of Figure 2, respectively. We name the two stellar components A and B, where A is 
the brighter component that lies closer to the galaxy (separated by ~0'/26), component G. 
The separation of the two stellar components derived from the i-band image is A6 = l'/14, 
and our relative astrometry is summarized in Table 2. The stellar components have similar 
colors of % — H ~ 1.6 and ~ 1.8 mag for A and B, respectively, while the galaxy G has a 
redder color of i — H ~ 2.4 mag. The photometry of the components is also summarized in 
Table 2. These astrometric and photometric properties are typical of known lensed quasar 
systems (e.g., Inada et al. 2006), and imply that both components A and B are lensed images 
of a quasar and component G is the lens galaxy. 



4. Spectroscopic Observations 

We carried out a spectroscopic observation of SDSS J1332+0347 using the Faint Object 
Camera And Spectrograph (FOCAS; Kashikawa et al. 2002) installed on the Subaru 8.2- 
meter telescope on 2003 June 20 with an exposure time of 900 sec in 0'.'8 seeing. The 
observation was conducted in 2x2 on-chip binning mode, using a O'.'6-width slit aligned 
along components A and B, with the grism 300B and the filter L600. This configuration 
provides a spectrum covering 3,900 A to 6,000 A with a spectral resolution of R~ 700, a 
spectral dispersion of 2.7 A pixel" 1 and a spatial pixel scale of 0'.'208. We used IRAF to 
analyze the spectra, extracting two traces corresponding to a blend of image A with the 



3 IRAF is the Image Reduction and Analysis Facility, a general purpose software system for the reduction 
and analysis of astronomical data. IRAF is written and supported by the IRAF programming group at 
the National Optical Astronomy Observatories (NO AO) in Tucson, Arizona. NO AO is operated by the 
Association of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the 
National Science Foundation. 
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galaxy G and image B, respectively. 

The two spectra are shown in Figure 3. We clearly see broad C III] emission lines 
redshifted to z — 1.4. On the red wings of the C III] emission lines, we also see an absorption 
doublet at 4686 A and 4728 A, which we interpret as Ca II H and K absorption lines 
redshifted to z — 0.191. These strong absorption lines significantly distort the emission 
lines. It would be helpful to have a cleaner emission line in order to determine the quasar 
redshifts more accurately, but, the C III] emission lines are the only strong emission lines 
in the FOCAS spectra. We estimated the quasar redshifts using three approaches; 1) using 
the whole line profiles (A Q b s ), 2) using only the line peaks (A pea k), and 3) using the whole 
line profiles after masking and then interpolating through the absorption lines (A mas k). The 
results, summarized in Table 3, confirm that the two quasars have the same redshifts. We 
adopt procedure 3) for our standard result of z — 1.445 as it is likely to be the most reliable. 
We also verified that the spectrum A+G has little contamination from quasar B by extracting 
a spectrum at the same distance from component B as component A+G but in the opposite 
direction. The spectrum, shown in Figure 3, has no significant emission, and we conclude 
that the C III] emission line in the spectrum A+G is due to component A. Furthermore, 
in the spectra of both components A+G and B, we can see Mg II AA2796,2803 doublet 
absorption lines and a Fe II A2600 absorption line by an intervening object at z = 1.119. 

The spectral flux ratio between components A+G and B, displayed in the bottom panel 
of Figure 3, indicates that combination of A+G is redder than component B. If we scale the 
spectrum B by the iJ-band flux ratio of 1.127 between components A and B, and subtract it 
from the spectrum A+G, then we obtain the residual spectrum shown in Figure 4. While the 
spectrum is noisy, especially around the C III] emission line and towards the blue wavelength 
region, the spectral energy distribution (SED) is consistent with a template spectrum of an 
elliptical galaxy (Kinney et al. 1996) at the redshift of the Ca II absorption feature, z = 0.191. 
Given the strength of the absorption, we conclude that the SEDs of components A and B 
are identical and that component G is an elliptical galaxy at z — 0.191. 

5. Lens Model 

We used the lensmodel package (Keeton 2001b) to fit a mass model to the observations 
of SDSS J1332+0347, assuming that the object is indeed a lensed quasar system. We assume 
the standard Singular Isothermal Ellipsoid (SIE) model characterized by eight parameters: 
the Einstein radius Re, ellipticity e, position angle 9 e , position of the lensing galaxy, the 
position and flux of the source quasar. Since the number of observable constraints is also 
eight (the positions of A, B, and G, and the fluxes of A and B), the number of degrees of 
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freedom is zero. As expected, we were able to find a model that perfectly reproduces the 
observables, x 2 ~ 0. The results of the model are presented in Table 4. We note that the 
best-fit values for the ellipticity and the position angle, e = 0.53±0.07 and 9 e = 25°±2°, are 
in good agreement with those observed for the galaxy G (derived using GALFIT in i-band 
image; see Figure 2), of e = 0.69 and 9 e = 21°, respectively. 

Einstein radii are related to the velocity dispersions of lensing galaxies, so we can 
use the Faber and Jackson relation (Faber & Jacson 1976) to estimate the luminosity of the 
lens galaxy. From the best-fit value = 0M65, we estimate that the apparent magnitude 
of the lens is % — 18.8, assuming the correlation of velocity dispersions and magnitudes of 
early-type galaxies derived by Bernardi et al. (2003). This is again in good agreement with 
the observed magnitude of the galaxy G, % = 18.6 (Table 2), implying that galaxy G is 
responsible for the most of the lens potential. 

6. Conclusions 

We have presented the extensive follow-up observations of SDSS J1332+0347 using the 
Subaru telescope and UH88 telescope. We found a bright galaxy (component G) between two 
stellar components (components A and B) with similar color in both the Subaru/Suprime- 
Cam i-band and UH88/QUIRC iJ-band images. Although the C III] emission lines of the 
two stellar components are distorted by the Ca II H and K absorption lines from the bright 
galaxy at z\ = 0.191, we can still confirm that the two quasars are at the same redshift, 
z s = 1.445. Subtraction of a scaled spectrum of component B from the component A+G 
leaves a residual that is consistent with the spectrum of an elliptical galaxy redshifted to 
z = 0.191. The observed luminosity, ellipticity and position angle of the bright galaxy 
are in good agreement with those expected from a standard lens model. All these results 
lead us to the conclusion that SDSS J1332+0347 is indeed a lensed quasar system, where 
a source quasar at z s = 1.445 is lensed by a bright lens galaxy at z t = 0.191 resulting 
in an image separation of A6 = l'/14. One peculiar characteristic of the lens is that the 
brighter component A is closer to the lens galaxy than the fainter component B. We have 
no difficulty reproducing the flux ratio, and note that such inverted flux ratios were also 
observed and easily modeled in the two-image lens HE1104-1805 (Wisotzki et al. 1993). 
This discovery adds another object to the complete sample of the SDSS lensed quasars that 
can be used to estimate the cosmological model. The brightness of the lens galaxy also makes 
it a good candidate for dynamical observations; higher-resolution imaging would allow us 
to determine more accurately the brightness of the components in this lens system, and 
spectroscopic observations with better signal-to-noise ratios would also allow us to study the 
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interstellar medium in the lens galaxy in detail. 
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Table 1. 



SDSS J1332+0347: SDSS PHOTOMETRY AND REDSHIFT 



u a g a r a i a z a Redshift b 

19.12±0.03 18.61±0.03 18.26±0.02 17.95±0.02 17.84±0.02 1.438±0.003 

a PSF magnitudes from the SDSS imaging data. 
b quasar redshift from the SDSS spectrum. 



Table 2. SDSS J1332+0347: ASTROMETRY AND PHOTOMETRY 



Object 


x[arcsec] a 


2/[arcsec] a 


i h 


H h 


A 


0.000 ±0.006 


0.000 ±0.006 


19.28 ±0.03 


17.72 ±0.17 


B 


-1.024 ±0.006 


-0.501 ±0.006 


19.66 ±0.02 


17.85 ±0.05 


G 


-0.228 ±0.013 


-0.119 ±0.013 


18.64 ±0.02 


16.28 ±0.04 



Note. - All the values are the outputs of the GALFIT procedure for 
SDSS J1332±0347. 

a The positive directions of x and y are West and North, respectively. 
The values are measured in the Subaru Suprime-cam i-band image. 

b The errors do not include the absolute calibration uncertainty. 



Table 3. SDSS J1332±0347: CENTRAL WAVELENGTHS OF C III] 



Object 


Aobs 


^obs 


^pcak 


%>eak 


-^mask 


^mask 


A+G 
B 


4653.49 
4663.26 


1.4380 
1.4431 


4661.11 
4666.08 


1.4420 
1.4446 


4666.23 
4666.25 


1.4447 
1.4447 
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Table 4. SDSS J1332+0347: LENS MODEL 



Parameter 


Value 


Einstein radius 


R E = 0"465 


Ellipticity 


e = 0.53 ±0.07 


Position angle 3. 


9 e = 25° ± 2° 


Total magnification 


that = 4.1 


Time delay- 


At = 6.7/i _1 day 



a Each position angle is measured East 
of North. 
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Fig. 1. — upper panel: The wide-field SDSS 2-band image of SDSS J1332+0347. The pixel 
scale is 0'/396 pixel -1 . In the image, North is up and East is left, lower panel: The SDSS 
spectrum of SDSS J1332+0347. 
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Fig. 2. — Follow-up images of SDSS J1332+0347. The optical i-band images taken with the 
Suprime-Cam on the Subaru telescope are shown in the left panels, while the near-infrared //-band 
images taken with the QUIRC on the UH88 telescope are displayed in the right panels. The original 
images, images after subtracting two PSF components, images after subtracting galaxy model, and 
images after subtracting both PSF and galaxy components are shown in the top, second, third, 
and bottom panels, respectively. In all panels, the box size is 20", and North is up and East is left. 
The pixel scales are Of! 202 pixel -1 (Suprime-Cam) and Of! 189 pixel -1 (QUIRC), respectively. 
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Fig. 3. — Spectra (R~700) of components A+G (thin black solid line) and B (thick black 
solid line) obtained with the Subaru/FOCAS in the upper panel. The spectrum of component 
A+G is shifted by constant to make it visible. A spectrum extracted at the same distance 
from component B as component A+G, but in the opposite direction, is also shown by the 
gray solid line. C III] emission lines (vertical solid line) in components A and B at z = 1.445, 
and Ca II H and K absorption lines (vertical dashed line) of galaxy G at z = 0.191 are 
detected. Fe II and Mg II absorption lines from an intervening system at z = 1.119 are also 
detected (vertical dot-dashed line). The Mg II absorption line doublet ratio of 1.0 is typical 
of such systems (e.g., Steidel et al. 1992). The spectral flux ratio between components A+G 
and B is shown in the bottom panel. 
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Fig. 4. — A spectrum of component G (gray solid line) obtained by subtracting a scaled 
(multiplied by 1.127) spectrum of component B from that of component A+G. The derived 
galaxy spectrum does not clearly show the Ca II H and K absorption lines, but is consistent 
with a template spectrum of an elliptical galaxy (Kinney et al. 1996) at z = 0.191 (black 
solid line). This result indicates that components A and B have identical SEDs and redshifts. 



